Clustering induced suppression of ferromagnetism in diluted magnets 
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Ferromagnetism in diluted magnets in the compensated regime p <C x is shown to be suppressed 
by the formation of impurity spin clusters. The majority bulk spin couplings are shown to be 
considerably weakened by the preferential accumulation of holes in spin clusters, resulting in low- 
energy magnon softening and enhanced low-temperature decay of magnetic order. A locally self- 
consistent magnon renormalization analysis of spin dynamics shows that although strong intra- 
cluster correlations tend to prolong global order, T c is still reduced compared to the ordered case. 



Diluted magnetic semiconductors (DMS) such as 
Gai-^Mn^As exhibit a novel carrier-induced ferromag- 
netism where doped mobile carriers mediate magnetic 
interactions between the S = 5/2 Mn ++ impurity 
spins with magnetization and Curie temperature 
T c increasing with carrier concentration^i^i^ While 
as-grown samples of Gai_ x Mn x As with high Mn concen- 
tration exhibit a large degree of compensation, where the 
hole density p is a small fraction of the Mn impurity con- 
centration x, recent progress in fabrication and annealing 
has yielded nearly compensation free samples, especially 
at low Mn concentrations. 

Strong thermal decay of magnetization observed in 
these diluted magnets has been characterized recently in 
terms of spin- wave properties*^ Magnetization studies of 
Gai-^Mn^As samples with different Mn content (thick- 
ness about 50nm, Mn content ranging from 2% to 6%) 
using SQUID (superconducting quantum interference de- 
vice) measurements, with linear extrapolation from a 0.3- 
0.4 T magnetic field to overcome anisotropy fields which 
result in spin-reorientation transitions^ have found the 
spontaneous magnetization to be well described by the 
Bloch form M (T) = A/ (1-BT 3 / 2 ), with a spin-wave pa- 
rameter B ~ 1 — 3x 10~ 3 K~ 3 / 2 which is about two orders 
of magnitude higher than for Fe and FeCo filmsi n i 13 The 
calculated values of spin-wave stiffness constant D were 
found to be of same order of magnitude as obtained from 
magnetic Kerr measurements using pump-probe setup of 
standing spin waves in Gai_ x Mn x As thin films^ 

Recent theoretical studies of magnon properties in the 
diluted ferromagnetic Kondo lattice model (FKLM) pro- 
vide quantitative understanding of this large magnitude 
of B in terms of a strong enhancement in density of low- 
energy magnetic excitations at higher dilution^ The re- 
sulting strong thermal decay of magnetization was found 
to be in good agreement with the Bloch form at low tem- 
perature, with B of same order of magnitude as obtained 
in the squid magnetization measurements. 

Generally, presence of competing antifcrromagnctic 
(AF) interactions is a source of magnon softening, due to 
the negative-energy contributions from AF bonds when 
spins are twisted away from the collinear ferromagnetic 
arrangement. Indeed, close to the limit p ~ x of half fill- 
ing, competing AF interactions become strong enough to 
destabilize the ferromagnetic state in the diluted FKLM. 



In the compensated regime p x, however, the relevant 
mechanism is quite different. 

In this paper we highlight a different aspect of com- 
peting interactions in diluted magnets — a competition 
between cluster and bulk spin couplings. In the com- 
pensated regime p <C x, the formation of spin clusters 
(due to impurity positional disorder) and accumulation 
of doped holes in these impurity-rich regions significantly 
deprives the majority bulk spins of holes. Because of the 
crucial compensation condition p <C x, even if cluster 
spins constitute a small fraction of the total, they can 
accumulate a majority of doped holes. While the cluster 
couplings are therefore dramatically enhanced, yielding 
high-energy cluster-localized magnon modes and strong 
intra-cluster ordering, the couplings between bulk spins, 
which are essentially responsible for long-range ferromag- 
netic order, are significantly weakened, resulting in low- 
energy magnon softening and enhanced thermal decay of 
magnetization. These effects of impurity clustering on 
spin couplings and dynamics have been qualitatively dis- 
cussed earlier within the impurity-band representation, 
applicable near the metal-insulator transition. — 

Our objective here is to quantitatively connect, within 
a more general and unified framework, the different rel- 
evant aspects — disorder, clustering, fermion polariza- 
tions, spin couplings, magnons, and spin dynamics. With 
respect to magnetic effects of impurity clustering in DMS, 
as distinguished from comparison with the ordered case, 
earlier studies have found softening of low-energy magnon 
modes, weakening (strengthening) of bulk (cluster) cou- 
plings, and formation of high-energy cluster-localized 
modeS ) 16 ' 17 i 18 formation of ferromagnetic droplets but no 
significant effects on T c in Monte Carlo calculations ; 1 ^ 20 
and reduction of T c in density functional^ and local spin- 
density approximation^ studies involving ab-initio eval- 
uation of exchange parameters and subsequent statistical 
analysis. 

The interplay between itinerant carriers in a partially 
filled band and localized impurity moments is conven- 
tionally studied within the diluted FKLM 

H = t ^2 a la a t+S,<r - ^J2 Sl - (Tl (1) 

involving itinerant electrons hopping on host sites i and 
a local exchange coupling between the localized impurity 
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FIG. 1: (Color online) Distribution of spin couplings (single 
configuration) at two doping concentrations, showing large 
spread in the strongly enhanced cluster couplings and dis- 
tinctly non-RKKY nature in the compensated regime (p -C 
x). 



spin S/ and itinerant electron spin <xj. Recently, magnon 
properties in the mixed spin-fermion model (1) were dis- 
cussed extensively with respect to variations in dilution 
x, exchange coupling J, impurity energy e^, and hole 
concentration p within an exact treatment of impurity 
positional disorder amd a non-perturbative treatment of 
exchange coupling^ The magnon propagator was stud- 
ied in the random phase approximation (RPA) , where the 
bubble diagrams, representing repeated interactions be- 
tween impurity spins mediated by the particle-hole bub- 
ble, provide the lowest-level spin-rotationally-symmctric 
treatment of transverse spin fluctuations where the Gold- 
stone mode is explicitly preserved. 

Due to spin-rotation symmetry of (1), the impurity 
spin dynamics is naturally described by an isotropic inter- 
action — Tin j JijSi-Sj between impurity spins, and at 
the RPA level, the (weakly-dynamical) carrier-mediated 
spin couplings are given by 



j ij {u) = L [x \ u))]ij 

in terms of the particle-hole bubble 
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evaluated by integrating out the fcrmions in the broken- 
symmetry state. Spin-fluctuation corrections to spin 
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FIG. 2: (Color online) Substantially weakened bulk spin cou- 
plings at separations > 2 compared to the ordered case, for 
a N = I0 3 system with J = 4, x = 1/8, and p = 1.35%, 
averaged over 20 configurations. 



couplings can be incorporated by including self-energy 
and vertex corrections in the particle-hole bubble within 
a spin-rotationally-symmetric scheme^ but are sup- 
pressed by the factor 1/5. The above isotropic spin 
interaction neglects magnetic anisotropy arising from 
spin-orbit interaction in the host semiconductor , 24 i 25 i 26 i 27 
and the recently observed strain-induced uniaxial 
anisotropy which is dependent on hole concentration and 
temperature^ 

For simplicity, we consider (1) on a simple-cubic lat- 
tice, with periodic boundary conditions. The N m mag- 
netic impurities are placed randomly on a fraction (/) 
of the N — L 3 host sites (i), with impurity concentra- 
tion x = N m /N . We consider positive nearest-neighbour 
hopping t, and set t = 1 as the unit of energy scale. 
For the fully polarized, collincar ferromagnetic ground 
state at T = 0, stability of which is confirmed from 
absence of negative-energy magnon modes, the fermion 
eigenvalues {Ei a } and wave functions {V'Zcr} are obtained 
by exact diagonalization of the N x N fermion Hamil- 
tonian with effective impurity potentials for the 

two fermion spins. Furthermore, we consider the sat- 
urated ferromagnetic state with a fully occupied spin-| 
band and hole doping only in the pushed-up spin- J, band 



(p = N^ oles /2N), so that only the first term in Eq. (4) 
contributes to the particle-hole bubble. For p ~ x or 
lower values of J, x, we do obtain negative-energy magnon 
modes, indicating noncollincar ferromagnetic ordering, as 
found in earlier studies^ As magnon energies are very 
low compared to the Stoner gap, the spin couplings are 
only weakly dynamical, and hence we have set u) = in 
Eq. (2). Although their nearly static nature is similar 
to that of RKKY interaction, the effective spin couplings 
are strongly non-perturbative in character due to the J, x 
dependence of the wavefunctions and eigenvalues in Eqs. 
(3) and (4). 
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FIG. 3: (Color online) Correlation between impurity-site 
fermion polarization (erf) and the local measure Jj = 
^2 Jij of spin couplings, obtained for 20 configurations 
of a N = 10 3 system with J = 4, x = 1/8, and p = 1.35%. 
Relative to the ordered case (bold dot) , significantly enhanced 
(reduced) fermion-spin densities and impurity-spin couplings 
correspond to cluster (bulk) spins. 



Figure 1 shows the distribution of impurity spin cou- 
plings Jij obtained from (2,4) as a function of impurity 
separation r = \rj — rj|, for two different hole concen- 
trations p/x w 1/10 and « 1/4. Relatively very strong 
spin couplings are generated between cluster spins having 
small separations 1 and y/2, which account for the high- 
energy cluster-localized magnon modest A most strik- 
ing feature is that for the same impurity separation, the 
coupling magnitudes are spread over a large range. The 
coupling between two impurity spins is thus not merely 
a function of their separation, but actually depends on 
the whole configuration, suggesting shades of a complex 
system^ Indeed, for same x,p,r, the bulk spin couplings 
in the disordered case are strongly weakened compared 
to the ordered case. This feature is absent in certain 
spin-only models where the couplings arc fixed functions 
of the separation, quite independent of disorder. 

Furthermore, while AF couplings are distinctly present 
for p/x m 1/4 and indeed the collinear ferromagnetic 
state is unstable, in the compensated regime p/x m 1/10, 
where the ferromagnetic state is quite stable, the cou- 
plings are seen to be dominantly ferromagnetic. The 
plot of configuration-averaged couplings (Fig. 2) con- 
firms that the couplings remain ferromagnetic and do 
not change sign. This distinctly non-RKKY behaviour of 
spin couplings in diluted magnets has been highlighted 
recently in a detailed study^ 

In order to examine the effect of impurity clustering 
on bulk spin couplings we make a comparison with an 
ordered diluted case. Fig. 2 shows the spin couplings 
for disordered (configuration-averaged over 20 configura- 
tions) and ordered cases for the same impurity concentra- 
tion and hole doping. Here the ordered case corresponds 
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FIG. 4: (Color online) The strong low-T decay of magnetiza- 
tion, with a subsequent marginally paramagnetic behaviour 
correspond to the enhanced low-energy magnon DOS and ten- 
dency of intra-cluster correlations to prolong order, shown for 
5 disorder configurations of a N = 12 3 system in comparison 
with the ordered case at exactly same dilution x = 1/8. 



to a superlattice arrangement of impurities with impurity 
spacing 2 and concentration x = 1/8. The bulk spin cou- 
plings involving average separations > 2 are clearly con- 
siderably weakened, which is directly responsible for the 
observed low-energy magnon softening with dilution^ 

Turning now to include fermions, and the main result 
of this paper, Figure 3 highlights the essential character 
of carrier-induced spin couplings in the presence of disor- 
der in terms of a direct correlation between local fermion 
polarization (erf) and the local measure Ji = yij-tj Jij 
of spins couplings. Fermion polarization values larger 
than the ordered-case value correspond to cluster sites 
which attract a disproportionately higher concentration 
of holes, yielding proportionately larger spin couplings. 
The remaining bulk sites are thus left with a hole deficit, 
resulting in proportionately weaker spin couplings. The 
compensation condition p<i plays a crucial role here 
because even if cluster spins constitute a small fraction of 
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the total, they can accumulate nearly all the doped holes. 
The three clearly separated cluster patches in Fig. 3 pre- 
sumably correspond to one-, two-, and three-dimensional 
clusters, with progressively higher degrees of derealiza- 
tion and impurity-level band broadening. 

The disproportionate accumulation of holes in 
impurity-rich cluster regions is simply due to preferen- 
tial filling of cluster impurity states by doped holes. 
Cluster impurity states undergo greater band broadening 
due to stronger effective hopping, and are thus pushed 
to the valence-band top. Impurity states correspond- 
ing to remaining bulk spins lying deeper in the valence 
band are left with hole deficit. With increasing J and 
more pronounced impurity-level character, this tendency 
of differential band broadening is enhanced, resulting in 
stronger suppression of fcrromagnctism, highlighting the 
non-perturbativc character of spin couplings in contrast 
to the RKKY picture 

To highlight the macroscopic effects of the disorder 
induced changes in bulk and cluster spin couplings, we 
compare the magnon density of states (DOS) and tem- 
perature dependence of lattice-averaged magnetization 
with the ordered case at the same dilution. Magnetiza- 
tion was obtained using a locally self-consistent magnon 
renormalization scheme,— which incorporates the spatial 
character of magnon states as well. Within the Tyablikov 
decoupling formalism, similar self-consistent schemes in- 
volving both global 3 ^ and locate self-consistency have 
been employed recently in the context of T c calculations 
for Heisenberg models with realistic spin couplings. 

Fig. 4 (upper panel) shows the characteristic disorder- 
induced features of sharply enhanced low-energy magnon 
DOS and high-energy structures due to cluster-localized 
modes, associated with the weak (bulk) and strong (clus- 
ter) couplings, respectively. The corresponding features 
in the magnetization plot (lower panel) are readily identi- 
fied. The rapid low-T decay follows from the low-energy 
magnon softening, and the subsequent marginally para- 
magnetic behaviour (more pronounced in some special 
configurations) is due to strong intra-cluster correlations; 
the overall behaviour is characteristic of two distinct 
(magnon) energy scales^ 

In Fig. 4 (lower panel) we have considered a bandwidth 
W = 12t = lOcV of the order of that for GaAs. For 
J = it = W/3 « 3.3cV as considered in our calculations 
and a nominal impurity concentration x = 10%, the value 



Jx = 0.33eV is close to the corresponding experimental 
value Jn Mn = J x 10% x n Ga = 150meV.A 3 xl0% x 
2 x 10 22 /cm 3 = 0.3cV with J = 150mcV.A 3 taken from 
Ref. [1]. Direct measurements yield J = 1.2 ± 0.2eV 
from core-level photoemission 3 ^ and J = 2.4±0.9eV from 
magneto-transport 

In conclusion, we have quantitatively studied the inter- 
play, within a non-pcrturbative and unified framework, 
of the different ingredients of ferromagnctism in diluted 
magnets — disorder, clustering, fermion polarizations, 
spin couplings, magnons, and spin dynamics. Due to the 
competition between impurities for a small fraction of 
doped holes in the compensated regime p <§; x, impurity 
disorder plays a highly non-pcrturbative role. The direct 
correlation between the broadly distributed fermion spin 
polarization values with the local measure of spin cou- 
plings highlights the essential character of carrier- induced 
spin couplings in the presence of disorder, clearly distin- 
guishing between cluster and bulk spins. Cluster cou- 
plings show a large spread for same impurity separation, 
and bulk couplings show distintly non-RKKY character. 

Preferential accumulation of holes in impurity clusters 
strengthens (weakens) the cluster (bulk) spin couplings, 
resulting in softening of low-energy magnons extended 
over bulk spins and formation of high-energy cluster- 
localized modes. The locally self-consistent magnon 
renormalization scheme is able to capture these disorder- 
induced features in the finite temperature spin dynam- 
ics. While the low-energy magnon softening is respon- 
sible for the observed strong low-T decay of magnetiza- 
tion (B w 1.6 x 10~ 3 K~ 3 / 2 in Fig. 4), strong intra- 
clustcr correlations tend to prolong global order through 
a marginally paramagnetic temperature regime. Inter- 
estingly, the somewhat slower decay due to cluster corre- 
lations brings the magnetization fall off closer to the T 3 / 2 
form over a broader T regime, as observed.— Overall, we 
find that ferromagnetism is considerably weakened com- 
pared to the ordered case. As increasing p/x enhances 
AF couplings and destabilizes the ferromagnetic state, 
reducing spin clustering appears crucial for enhancing 
ferromagnetism in diluted magnetic semiconductors, and 
indeed could be a contributing factor in the annealing 
processiii A two-component picture of bulk and cluster 
spins, with a relatively low carrier density in the bulk 
could be useful in interpreting the relatively high resis- 
tivity seen even in the most metallic samples 



* Electronic address: avinas@iitk.ac.in 

1 F. Matsukura, H. Ohno, A. Shen, and Y. Sugawara, Phys. 
Rev. B 57, R2037 (1998). 

2 T. Dietl, H. Ohno, F. Matsukura, J. Cibert, and D. Fer- 
rand, Science 287, 1019 (2000). 

3 H. Ohno and F. Matsukura, Solid State Commun. 117, 
179 (2003). 

4 C. Timm, J. Phys.: Condens. Matter 15, R1865 (2003). 

5 I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 



76, 323 (2004). 

6 H. Boukari, P. Kossacki, M. Bertolini, D. Ferrand, J. Cib- 
ert, S. Tatarenko, A. Wasiela, J. A. Gaj, and T. Dietl, 
Phys. Rev. Lett. 88, 207204 (2002). 

7 K. C. Ku, S. J. Potashnik, R. F. Wang, S. H. Chun, P. 
Schiffer, N. Samarth, M. J. Seong, A. Mascarenhas, E. 
Johnston-Halperin, R. C. Myers, A. C. Gossard, and D. 
D. Awschalom, Appl. Phys. Lett. 82, 2302 (2003). 

8 K. W. Edmonds, P. Boguslawski, K. Y. Wang, R. P. Cam- 



5 



pion, S. N. Novikov, N. R. S. Farley, B. L. Gallagher, C. 
T. Foxon, M. Sawicki, T. Dietl, M. B. Nardelli, and J. 
Bernholc, Phys. Rev. Lett. 92, 037201 (2004). 
A. M. Nazmul, S. Kobayashi, S. Sugahara, and M. Tanaka, 
Jpn. J. Appl. Phys. 43, L233 (2004). 

10 X. Liu, W. L. Lim, M. Dobrowolska, J.K. Furdyna, and T. 
Wojtowicz, Phys. Rev. B 71, 035307 (2005). 

11 M. Sperl, J. Sadowski, R. Gareev, W. Wegscheider, D. 
Weiss, and G. Bayreuther, to be published. 

12 M. Sawicki, K.-Y. Wang, K. W. Edmonds, R. P. Campion, 
C. R. Staddon, N. R. S. Farley, C. T. Foxon, E. Papis, E. 
Kamihska, A. Piotrowska, T. Dietl, and B. L. Gallagher, 
Phys. Rev. B 71, 121302(R) (2005); M. Sawicki, J. Magn. 
Magn. Mater. 300, 1 (2006). 

13 W. Kipferl, M. Sperl, T. Hagler, R. Meier, and G. 
Bayreuther, J. Appl. Phys. 97, 10B313 (2005); M. Sperl, 
W. Kipferl, M. Dumm, and G. Bayreuther, ibid 99, 08J703 
(2006). 

14 D. M Wang, Y. H. Ren, R. Merlin, K. Dziatkowski, X. 
Liu, J. K. Furdyna, and M. Grimsditch, to appear in ICPS 
proceedings (2006). 

15 A. Singh, S. K. Das, A. Sharma, and W. Nolting, 
|cond-mat/0607633| (2006) . 

16 M. Berciu and R. N. Bhatt, Phys. Rev. Lett. 87, 107203 
(2001); Phys. Rev. B 66, 085207 (2002). 

17 A. Singh, A. Datta, S. K. Das, and V. A. Singh, Phys. Rev. 
B 68, 235208 (2003)^ 



A. Singh, |cond-mat /030 7009 (2003). 

M. Mayr, G. Alvarez, and E. Dagotto, Phys. Rev. B 65, 
241202(R) (2002). 

D. J. Prior, Jr. and S. Das Sarma, Phys. Rev. B 73, 165203 
(2006). 



21 L. M. Sandratskii, P. Bruno, and S. Mirbt, Phys. Rev. B 
71, 045210 (2005). 

22 J. L. Xu, M. van Schilfgaarde, and G. D. Samolyuk, Phys. 
Rev. Lett. 94, 097201 (2005). 

23 A. Singh, Phys. Rev. B 74, 224437 (2006). 

24 J. Schliemann, J. Konig, and A. H. MacDonald, Phys. Rev. 
B 64, 165201 (2001). 

25 G. Zarand and B. Janko, Phys. Rev. Lett. 89 047201 

(2002) . 

26 L. Brey and G. Gomez-Santos, Phys. Rev. B 68, 115206 

(2003) . 

27 M. Berciu and R. N. Bhatt, Phys. Rev. B 69, 045202 

(2004) . 

28 J. Schliemann and A. H. MacDonald, Phys. Rev. Lett. 88, 
137201 (2002). 

29 N. Goldenfeld and L. Kadanoff, Science, 284, 87 (1999); 
E. Dagatto, Science, 309, 257 (2005). 



R. Bouzerar, 



Bouzerar, and T. Ziman, 



|cond-mat/06 07640 (2006). 

31 G. Alvarez, M. Mayr, and E. Dagotto, Phys. Rev. Lett. 
89, 277202 (2002). 

32 S. K. Das and A. Singh, |concTmat/0506523l (2005). 

33 S. Hilbert and W. Nolting, Phys. Rev. B 71, 113204 (2005). 

34 R. Bouzerar, G. Bouzerar, and T. Ziman, Phys. Bev. B 
73, 024411 (2006). 

35 J. Okabayashi, A. Kimura, O. Rader, T. Mizokawa, A. 
Fujimori, T. Hayashi and M. Tanaka, Phys. Rev. B 58, 
R4211 (1998). 

36 T. Omiya, F. Matsukura, T. Dietl, Y. Ohno, T. Sakon, M. 
Motokawa, and H. Ohno, Physica E (Amsterdam) 7, 976 
(2000). 



